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Abstract We consider N nearest neighbor random walks on the positive integers with a drift 
towards the origin. When one walk reaches the origin, it jumps to the position of one of the other 
— 1 walks, chosen uniformly at random. We show that this particle system is ergodic, and 
establish some exponential moments of the rightmost position, under the stationary measure. 
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1 Introduction 

There is recent interest in approximating the limiting law of irreducible Markov processes condi- 
tioned not to hit some (forbidden) state [HI [15]. This limiting law is not guaranteed to exist, but 
when it does it is called a quasi-stationary distribution (QSD). For QSD in the context of Birth 
and Death chains, we refer to |9j, and the situation treated there is one in which there is a one 
parameter family of QSD. 

QSD are neither well understood, nor easily amenable to simulation. One proposal made by 
Burdzy, Holyst, Ingerman, and March (in a particular setting) is to consider independent 
Markov processes except that when one reaches the forbidden state, it jumps to the state of one 
of the other processes, chosen uniformly at random. The natural conjecture is that the empirical 
measure, under the stationary measure, converges to the QSD as the number of processes goes 
to infinity. It is also natural to conjecture that the selected QSD is the minimal., in terms of average 
time needed to reach the forbidden state. 

In this note, we consider random walks on N, in continuous time, with a drift towards the 
origin. When one random walk reaches the origin, it jumps instantly to the position of one of the 
other N — 1 walks, choosen uniformly at random. We call Fleming- Viot the interacting random 
walks just described. Indeed, this dynamics has a genetic interpretation. 



• Positions of the walks are evolving genetic traits of A^ individuals. 

• The forbidden state (here 0) is a lethal trait [the selection mechanism). 

• At the moment an individual dies, another one, chosen uniformly at random, branches {the 
branching mechanism). This keeps the population size constant. 
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We establish a Foster's criteria, which gives ergodicity, as well as a control of small exponential 
moments of the rightmost walk. To state our main result, let denote the position of the N 
interacting walks at time T, and let £'[-|^o = i] denote average with respect of the law of the 
process with initial condition ^. 

Theorem 1.1 There are positive constants K,a, k,5q, A,ci,C2,c^ such that for G N, time T = 
j41og(A^), any 5 < 60, and ^ G N^, we have 

E [exp (,5 max(eT)) jCo = (\- exp {5 max(e)) < - ci lmax(c)>Kiog(7V)e^°'^''^^^ 

j_ „ U „-KT„<5max(0 I „ „i5alog(Ar) ^ ' ' 

+ C2-"max(5)>i^log(Af)e ^ ^''^+036 '. 

As a consequence, for N large enough there is a unique invariant measure Xn for Fleming-Viot. 
Integrating over X^, there are /3, C > such that for any N, and 5 < 60 

y"exp(5max(0)dA7v(0 < Cexp (5/31og(A^)). (1.2) 

This first elementary step is an important ingredient in the proof of the conjecture we alluded 
to above [13] • Also, it might be of independent interest in view of recent deep and comprehensive 
studies on the rightmost position in branching random walks [U [U [101 HI EJ [H O [12] . This selection 
of recent works is far from being exhaustive, but already shows the vitality of this issue. 

The rest of paper is organized as follows. In Section [2] we define the model, and recall well- 
known large deviations estimates. In Section [3l we explain how to divide walks into groups with 
little correlations over a well chosen time period. In Section HI we estimate the probability that the 
maximum displacement does not decrease. Finally, in Section [Sj we establish Foster's criteria. 



2 Model and Preliminaries 

Here, we deal with continuous-time nearest neighbor random walks on N, with rate p to jump right, 
and rate q = l— p>p to jump left. The drift is —v with v = q — p > 0. A single walk makes Nt 
jumps in the time period [0,t], and its increments are denoted Xi, . . . ,Xn, with = —v, and 

Xi = Xi + V denotes the centered variable. Note that 

T 

P(^(Xi + i;) > xT) < exp(-r/(x)), (2.1) 
1=1 

with 

I{x) = sup {Xx - A(A)} with A(A) = log(pe^ + qe~^) + Xv. (2.2) 

A>0 

Due to the nearest neighbor jumps of our walk, we have 

I{v + 1) = log — ) X > V + 1, I{x) = 00. (2.3) 

We define also x 1— )■ I{x) = 1 — exp(— /(x)), which is discontinuous at v + 1 with 

i{v) < i{v + l) = q and for x > v + 1, i{x) = 1. (2.4) 
Note that if Nt is Poisson of mean T, then 

P{ J2 i^i + ^)^ ^^t) < exp(-r/(x)). (2.5) 

i<NT 
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On Poisson tails. We need two rough tail estimates on the Poisson clocks. Both are obvious 
and well-known. Assume T are positive. 



P(A^T>er + x) <exp(-T-x), 



(2.6) 



and 



Pi^Nr <-T-x)< exp(-(l - 2/e)T - x), 



(2.7) 



Both are obtained readily by Chebychev's inequality. Indeed, we obtain (j2.6p from 

P{Nt >eT + x)< e-^^-^E[e^^] = exp(-r - x), 



(2.8) 



and we obtain (12. 7p from 



P(^Nt <-T-x)< e^''-^E[e-^^] = exp(-(l - 2/e)T - x)- 



(2.9) 



e 



3 Independence 

On the multitype branching of [6j . A key idea introduced in [6J is to embed the Fleming- Viot 
process into a multitype branching process whose space displacements and branching mechanism are 
independent, and which is attractive. We refer to Section 3 of [6j for a description of the multitype 
branching process, and recall here its main features. Assume that we start with interacting 
random walks. This defines types with which we associate N independent exponential clocks of 
intensity q with marks. The time realizations of the clock of type i have marks in the set of labels 
{1, . . . , A^}\{i}, and each mark is chosen uniformly at random from the A^ — 1 symbols. When 
clock i rings, and when its mark is j, each walk of type j branches into two children: one of type i 
and one of type j. The two children behave as independent random walks starting at the position 
of their parent. If B^^ denotes the population of individuals alive at time T, and |Bt| denotes its 
cardinal, it is easy to see the equality £'[|D'r|] = |Do| exp(qT). For an individual v alive at time T, 
we denote by 1 1— >■ S^it) its trajectory for t G [0,T]. 

Independent groups of walks. A drawback of the multitype branching process is an expo- 
nentially growing population. Since, we use a time of order log(A^), we cannot use here such an 
embedding. Even though in the Fleming- Viot process, all particles interact which each other, a sim- 
ple observation is that as long as a particle has not touched the origin its trajectory is independent 
from the other ones, even though this trajectory might influence others. 

To create some independence between walks, we decompose the interacting walks in two sets 
at time 0. We first fix a time T and a length L to be chosen later. 

• The blacks, whose initial position is below L. 

• The reds, whose initial position is above L. 

Then, color changes as follows: if a black walk jumps on a red walk, it becomes instantly red. We 
interpret this jump as a red binary branching. Now, red walks are not independent from black walks 
because they might touch the origin before time T, and jump onto a black position. However, if 
vT <C we expect this to be rare. To obtain independence, we add another color to our description: 
each red walk is coupled with a green walk which behaves identically in terms of move or branching 
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but with green children, except that when a green walk reaches the origin it continues its drifted 
motion on Z (without selection mechanism). Thus, green walks behaves like independent random 
walks with branching at the times a black particle hit zero and chooses the label of a green walk. 
If Rq is the first time one of the red walks touches the origin, we have that at time T, on the event 
{Rq > T}, red and green positions are identical. The point of introducing green walks is that 
their branching times is independent of their positions. We denote with D^, Dj,, the respective 
number of red, green and black walks at time T. Also Dt = DJ^ U D?p U = {!,..., A^}, and we 
still denote hy S^it), the trajectory of f G Dt- 

When embedding a group of walks into a branching process, we denote with B^,D|,,B|, the 
respective number of red, green and black individuals in the mutlitype branching processes. 

The key idea here is to work on a time of order log(A^), to control the black walks by a multitype 
branching process, but to let the red walks (or rather the green walks) grow as in Fleming- Viot 
with a population bounded by A^, and with branching due to independent black walks. 

On the choice of time T and length L. We choose T large enough so that q + log{N)/T < 1. 
We actually need a little more. 

. = min (1 - 2/e. /(|) - l^iM, 1 _ , _ l^lM) > 0. (3.1) 

Once T = Alog{N) satisfies dXU, we set L = eT. 

4 When things go wrong 

We wish to estimate the probability of the event where the maximum displacement does not de- 
crease. We thus define a had set B(T,L) as containing the following events: 

• One red walk reaches the origin before time T (i.e. {Rq < T}). 

• One black walk travels a distance L upwards in a period [0, T]. 

• The maximum displacement of a green walk in a time T is above —j^T. 

Thus, on the complement on B{T,L), green and red are identical, and 

max Sy(T) — max 5^,(0) < max (Sv(T) — 5^,(0) ) = max (Sy(T) — 5„(0) ) < T, 

vGDt vGDo v&Dt V / fGDf, V / 2e 



which implies that if M{T) = max^jg/jj, Sy{T) 

E [Ib.(t,l) exp (<5(M(T) - M(0))) |C(0) = e] < exp ( - ^T) . (4.1) 
We estimate next the probability of each event making up B{T,L), with the following outcome. 
Lemma 4.1 For any ^ S N^, we have, with k as in liS.l]) . 

P{BiT,L)\m = < 4exp(-/^r). (4.2) 
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4.1 A red walk does reach 



Recall that L = eT. We embed the Fleming- Viot into a branching multitype, while keeping the 
red coloring. We need to estimate the probability that one red displacement gets below L units in 
a time period [0,T]. Note that to realize {Rq < T}, there is v G such that the number of its 
time jumps Nt must be larger than L, and this is what we use. 

p{Ro < T\m = <ewt\] X p{3t <T, Y.Xi< -^^) 

i<Nt (4.3) 
<E[\BJr\] X P{Nt > eT) < Ne'^^e''^ < e~^^. 

4.2 A large black displacement 

Recall that at the time a black reaches 0, and jumps on a red walk, it ceases to be black to become 
red. We bound here the black walks with a multitype branching, assuming that blacks do only 
jump on blacks, with the effect that we are overestimating the black population. The estimates 
are similar to these of Section [4.11 We use that to make L steps right, a black walk must make L 
time-marks {Nt > L), and this event is estimated in (j4.3p . 

4.3 Green's maximum too high 

The key point is that the green branching times are independent of positions of the green. They 
depend only on the history of black walks. Also, the population of green walks is bounded by 
N. Thus, it is crucial here not to use the multitype branching of [6j: we estimate the probability 
that {max^g^s {Sv{T) — ^^.(O)) > —vT/{2e)}. Define Mt{i) as the number of green walks whose 
displacement during time period [0, T] is larger than 7. Then, 

EWT{i)\m = i] = E[\D?,\\m = e] X p( ^ X, > 7). (4.4) 

i<NT 

The reason is the independence of the branching times and displacements of the green walks. For 
V G Z?^, assume there are v branchings before time T, say at times Ti, . . . ,Tj^, and we have (for 

xf''' i.i.d. independent from {Ti,i € N}) 

S,{T)-SM= E + E ^"t^- (4.5) 

As one conditions first on the black history up to time T, one fixes the times Ti,...,Tj^, and 
obtain that ^[0, Ti] + • • • +N\Ty,T] sums up to a Poisson variable X[0, T] of intensity T, and most 
importantly 

E #^ + ---+ E E (4-6) 

ieAf[0,Ti] i&N[T^,T] i&N[0,T] 
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where the {Xi, i E N} are i.i.d. increments independent of iV[0,r]. We obtain, with k defined in 



P{ max(5,(r) - 5,(0)) > - -T\m = C) < E 



A/'H-f )|e(o) = c 



<NP{ X, > vNt - —T) 

i<NT 



<^(P{ Y.^i> I^t) + P{Nt < -T 



i<NT 

<iV(^exp ( - TlC^)) +exp ( - (1 - 2/e)r)) < 2exp(-Kr). 

5 Foster's criteria 

We start with an estimate on the tail, and of the exponential moments. 
5.1 On exponential moments 

We deal here with the multitype branching process. Recall that Bq = {1, . . . , N}, and let 5(0) 
{5„(0),?; GBo}. 



Lemma 5.1 For any T satisfying h3. and any X > 

P( max (5,(r) - 5,(0)) > eT + x|5(0) = i) < exp(-x). (5.1) 

Proof. Since the branching mechanism is independent of positions 

Nt 

P{ max (5,(r) - 5,(0)) > eT + x|5(0) = C) < E[\Bt\] x P( V X, > eT + x) • (5.2) 

i=l 

Now, if Xi denotes the centered variable, note that since the walk is nearest neighbor 

Nt 

> iv + l)NT) = 0. 

i=l 



Now, 



Nt Nt 
P{J2Xi>eT + x) =P{Y.X,>vNT + eT + x) 

i=l i=l 

Nt 

<P{ J2x,>{v + 1)Nt) + P{Nt >eT + x) 

i=l 

<0 + P{Nt >eT + x) 



(5.3) 



Now, Nj- is a Poisson variable of mean T, the standard estimate ()2.6p leads to 

P{Nt> eT + x) <ew{-T-x)- (5-4) 



Also, we have S[|Br|] < iVexp(gr), and ([O]) and the choice of T m yield 

P( max {S^{T) - 5^(0)) > eT + x|5(0) = i) < Ne'^^e'^'^ < e'^. 



(5.5) 



We can state our main estimate. 



Lemma 5.2 Assume that (1 — q)T > log{N), and 5 <1. Then, we have 



E 



exp f,5(niaxCi(r) - niaxC,(0))) ^(0) = ^ 



< 



1 



JeT 



1-5 



Proof. For any random variable X, we have 



E 



jx 



-l + / 5e^"P(X > u)du 
'o 



pel poo 

<1 + / de^'^du + / 5e'^"P(X > u)du 

Jo JeT 



<e'^"^ (1+1 Se'^'PiX >u + eT)du 



Now, using the tail estimate (|5.1|) . we have 



exp I 6 max (5^(T) — ^^((O)) 



E 



We now use the following bound to conclude 



JeT 
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(5.6) 



(5.7) 



(5.8) 



E 



exp ( 5(maxei(r) - maxCi(O)) ) 1^(0) = ^ <E exp ( (5( max 5„(T) - max5^(0)) ) |5(0) = ^ 

exp f 5 max (5„(T) - SJO))] \S(0) = ^ 



(5.9) 



5.2 Proof of Theorem [TT] 

We recall the general strategy of the proof of Proposition 1.2 of [6] (The Foster criteria). We have 
a bad set B{T,L) (which depends on T and L) which contains the cases where the maximum 
increases over a period [0, T], or when black walks win over or influence red ones. First, there is a 
set K on which we do not expect the maximum to decrease, with 



K = |max(5„(0)) < 3l} . 
Then, there is a good set where the maximum decreases: 



(5.10) 



K^r\B^{T,L) CG = {max(SJT)) - max(5„(0)) < -— tI . (5.11) 

[veDT v&Dt 2e J 
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Now, set Mf = max5t,(t). When is the initial configuration, and when we work with ^ E K'^, we 
have using Cauchy-Schwarz 



E 



<%i^ce'5*^«(P(i?(r,L)|C(0) = [exp(25(Mr - Mo))|C(0) = e] 
^gi<-ce""^" I 1 - e 



1/2 



1 



(5.12) 



^25eT 



We know from Lemma 15.21 that for 5 < 1/2 

E [exp(25(MT - Mo)) 1^(0) = i] < 
Note also that on the set K, if we use Lemma 15.21 



Thus, adding (|5.12p and (I5.13p . we obtain 



< Iifexp(3(5L + (5er). 



(5.13) 



E 



(5.14) 

We use now Lemma l4.ll and choose 5 small enough so that k > 4Je with the result 



E 



VI -2(5 



(5.15) 



Inequality (j5.15p is a Foster's criteria (see |13|, Theorems 8.6 and 8.13]). This implies the first part 
of Theorem 11.11 

Now, as we integrate (j5.15p with respect to the invariant measure, the left hand side of (|5.14p 
vanishes, and we obtain 

1 _ e-^^^/^) [ e^^^^Ux^iO < exp{36L + 5eT) + exp(-^r) / e^^(«)dA^(e). (5.16) 

With A large enough so that (j3.ip holds with T = yllog(A^) and L = eT, the second part of 
Theorem 11.11 follows at once. 

Acknowledgements. We would like to thank Djalil Chafai for valuable discussions. 
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